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In this paper we study the nonlinear current response of (7(1) spin liquids with large spinon Fermi 
surfaces under the perturbation of a time-dependent ac electric field E(t) within the framework 
of an effective U{1) gauge theory. In particular, the third-order nonlinear current response to ac 
electric fields is derived. We show that as in the case of linear current response, an in-gap power-law 
a;’’) response is found for the nonlinear current at low frequency. The nonlinear susceptibility 
may also induce through process of third harmonic generation propagating EM wave with frequency 
Suj inside the spin liquids. 


A Quantum spin liquid state is a Mott insulator with¬ 
out any long-range magnetic order down to the lowest 
temperatures [l|, For a long time, this exotic state 
exists only in the theorists’ models. Recent experimen¬ 
tal progress on the organic compound /v-(ET)^^(CN)3 
[j-Q and a few other candidate materials have 

changed the situation. In the case of k-(ET)2Cu 2(CN)3, 
the system is described by a nearly isotropic spin-1/2 
Hubbard model on the triangular lattice, and it is be¬ 
lieved that the insulator state is near the Mott transition, 
where the charge excitations are gapped and the spin-1/2 
spinons form a Eermi surface. The dynamics of the spin 
liquid state is described by an effective C/(l) gauge the¬ 
ory [HJ 1 ^ . Several experimental results have confirmed 
the Fermi-liquid-like behavior like the linear temperature 
dependence of the specihc heat, the Pauli-like spin sus¬ 
ceptibility, the nearly unity Wilson ratio, etc.. These ex¬ 
periments suggest a deep relationship between the unfa¬ 
miliar spin liquid state and the well-studied Fermi liquids 


In the continuum limit, all existing theoretical studies 
suggested that the spin liquids can be described by an 
effective 17 ( 1 ) gauge theory Lagrangian: 

: V’c'') (^0 i : 

C = il)l {idt - ao- Ao)ipcr - ■;r^|(V - fa - 

ZTTIq 

^a = Yg 4 (S) +Cf\ao) + Cf{ao,: 

Here 'i/J(^/cr) are spin-cr fermion creation (annihilation) 
operators, (Ao,A) is the external applied electromag¬ 
netic field, and g and Va are the coupling constant and 
effective velocity of the internal gauge field, respectively. 

is the effective (quadratic) action for the longitudinal 
gauge field oq. The dynamics of the internal gauge field 
arises from integrating out the charge degree of the sys- 
tem through for example, the slave-rotor approach 13 1 
or from an effective Landau Fermi-liquid-like approach 
[l 5 | . typically to second order in = d^a,, — d,ja^. In 
this paper we keep also the d**' order terms in Cg, 

£^■*^(00, a) ^ 6i(e-e)^-|-&2(b-b)^-|-63(e-e)(b-b)-|-64(e-b)^, 


( 1 ) 


where e = —dta— Voq and b = V x a are the internal 
electric and magnetic fields associated with the internal 
gauge field a^, respectively. These terms arises naturally 
in a lattice gauge theory model [I^. We have assumed 
that the spin liquid system respects time-reversal, inver¬ 
sion and rotational symmetry in writing down . The 
dynamics of internal gauge field are all expressed in terms 
of the field strengths e and b because of gauge invariance 
associated with internal gauge field in the insulator state 


Previous studies in the electrodynamic properties of 
spin liquid states mainly focused on the linear response 
of the system, i.e. the conductivity tensor and dielectric 
function. In this paper we study the nonlinear response 
of C/(I) spin liquids with large Fermi surface to external 
EM fields. We shall focus on the q —>■ 0 , w ^ 0 responses 
in this paper and shall derive the third-order electric con¬ 
ductivity for the spin liquid state. In particular, 

we are interested at examining whether a non-linear I—V 
relation I ^ 17 ^+^ (/3 > 0 ) may exist in spin liquid states. 
This is driven by the predicted existence of power-law 
linear conductivity aA')(uj) ~ in 17 ( 1 ) spin liquids. 
The power-law conductivity leads to nonlinear I—V re¬ 
lation in many one-dimensional Luttinger liquid systems 


Our approach consists of two steps: hrst we derive an 
equation of motion for the fermion (spinon) current un¬ 
der the perturbation of a time-dependent but spatially 
uniform electric field E(t) from which the nonlinear re¬ 
sponses of the spinons to the total electric field can be 
computed. The equation of motion is supplemented by 
an equation determining the internal electric field e and 
the two equations of motion are solved self-consistently 
to determine the nonlinear responses of the spin liquids 
to external electric fields. 


Equation of motion for fermion (spinon) current 

We first consider the equation of motion for the fermion 
current under a time-dependent but spatially uniform 
electric field. In this case it is most convenient to perform 








2 


a transformation to a non-inertial frame which moves 
with the center of the mass of the spin liquid system 1121. 


In this coordinate system, the time-dependent center-of- 
mass position vector R(t) = l/N'^-{ri{t)) {N = total 
number of fermions) satisfies the equation of motion: 

= j d'^rn(r,t)£:(r,t) - ^ f d^r n{r,t)U{r) 


+ 


1 J d‘^rjir,t)xBir,t), 


( 2 ) 


where n(r, t) = <5(r — r^)) is the time dependent par¬ 

ticle density distribution function. And 




nir,t) 

Too 


(A(r,t) +a{r,t)) 


is the total current density where jo(r,t) = 

Vi) + (5(r — fi) 2 ^]) canonical current density, p = 

canonical momentum, f = E -|- e and S = B -|- b are the 
total electric and magnetic fields (E and B are the exter¬ 
nal electric and magnetic fields, respectively), and U{r) 
is a background disordered scalar potential the fermions 
see. The meaning of this equation is clear: the first term 
on the right hand side is just the electric force, the second 
term comes from disordered potential and the last term 
corresponds to the Lorentz force contribution. 

The disordered potential U{r) = Xu{r) satisfies: 

(w(r)) = 0, (u(r)M(r')) = \u\'^S{r - r'), (3) 

and can be integrated out to second term in [13: 

- - W),t -1'), (4) 

where Xd(i" — — t') is the density-density response 

function for the system in the absence of impurities. 

The equation can be simplified further if the ap¬ 
plied electric field E(t) ^ Eoe“*‘^* is weak and it’s fre¬ 
quency uj is small. In this limit R(t) is slowly vary¬ 
ing and we may approximate R(t) — R(lO R(t)(t — 
t') + (I/2)R(t)(t — t')"^ and neglect higher order time- 
derivative terms. Then for a system with inversion sym¬ 
metry Xd(—q,w) = Xd(q,w), the right hand side of 
Eq. dH) is related to the density-density response function 
through [13: 

J dt' VR(4)Xd(R(t) - - t') 

= ^ *q ("l - • R(t)^^ Xd(q, w')L.=q.R(t). (5) 

q ^ ^ 

Putting together Eqs. (ED, O and ©, we obtain the 
equation of motion: 

m.R(t) = ^ J d^rn{r,t)£{r,t) ( 6 ) 


-V qimxd q, q • R.W) + 


N 


d''rj(r,t) X 


where the effective mass is given by: 


.A U x y C/ , J ^ I / pr\ 


where d = system dimension. Notice that the second 
order term R(t)(t — t')^ only leads to a mass renormal¬ 
ization that is proportional to the strength of disorder. 

As ImX(i(q,w) ^ 0 only at \uj/qvF\ < 1, the leading 
order contributions to the imaginary part of the density- 
density response function is given by: 

ImXd(q,w) ~ ai(g)—+ 02(9) ( — ) H-, ( 8 ) 

qvF \qvFJ 

where even order terms are zero, ai{q) and 02 ( 9 ) are 
w-independent parameters. In 2D, we obtain 


ai{q) = iV(0) 

a2(<?) = N{d) 


2 \2kF 


+ 0{{q/2kFf) 


2 d\2kF) 


+ 0{{q/2kFf) 


Thus we obtain 

J dt' VR(,)Xd(R(t)-R(t'),t-0 

= ——R(t) -|- 7TO^R(t)R(t)^, (9) 

T 

in the small w limit where r and 7 are related to ai{q) 
and 02 ( 9 ) by 

— = / dqdil q'^ai{q) cos^ 9, 

T nkF J 

7 = — ^ J 9'^“2((z) cos"^ 0. (10) 

Therefore in the absence of external magnetic field, the 
equation of motion for the charge current j(t) = noeR(t) 
in the weak field limit is: 

TTiiit) = noe^E(t) — — j(t) -I- noe^ J d'^r [no + Sn{r, t)]e(r, t) 

+ (^^) 2 jW[jW -jW] + e j d‘’*rj(r,t) x b(r,t), (11) 

where the particle density has been decomposed into a 
constant part no and a varying part 6n(r, t). In the q —>■ 0 
limit, 5n(r,t),b(r,t) —>• 0 and the equation reproduces 
the usual Drude conductivity if the non-linear term 
and the internal electric field e = — 9ta — Voq are ne¬ 
glected. 
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Third order nonlinear current response 

To solve the above equation we need to determine the 
internal gauge field a. The internal gauge field a can 
be determined by the Lagrangian equation of motion 
5Ces/5aL = 0 , which leads to the equation of motion (in 
the Coulomb gauge and in the limit q —?> 0); 


com pos ition rule in the framework of t/(l) gauge the¬ 
ory [IJ, ll5j, which leads to the power-law optical con¬ 
ductivity in the Mott gap. The nonlinear terms are 
results of non-linearity in the quasi-particle scattering 
(Imyd) and the presence of In the small-w limit, 

^ Lo/{ij3), and the nonlinear response reduces to 
a simple form: 




( 12 ) 

Eqs. (fTT|) and (fT^ together determine the current 
response under the time-dependent electric field E(t), 
which contain nonlinear effects up to third order in E. 
To solve the equations, we first Fourier transform Eq. 
(fT^ in the limit q — 0 , 


—e(w) -I- ^ e(wi)[e(wm) • e(w„)] = 

^ (Imn) 


E(c^0(E(w™)-E(w„))(.16) 

{Imn) 

Therefore our theory predicts that in the low frequency 
limit, there is a third-order correction to the current, i.e., 

j(a;) = cr(i)(a;)E(a;) 

+ (w;-|-Wm-|-u;„)E(aj;) (E(a;m) ■ E(a;„)) (17) 

{Imn) 

with proportional to ct^^^(w). In particular, we 

do not find any nonlinear I — V relation. 


wliere E(lmn) = -UJl-UJm- UJn). 

Therefore up to third order, the internal electric field 
e(a;) is related to the total current j(a;) by 


e(w) 


iuj ■’ ^ UJiUJrn^ri 


where b = Abig^v^. 

Lastly, we insert this result into Eq. (ED, obtaining 


/ \ • / \ 9 I / \ 9 9 1 f ^• / \ 

w(-ia;)j(a;) = noe E(a;) -f noe gv^- -j(a;) 

lUJ T 


+ 


7771'^ 

(noe)2 


E 

{Imn) 


1 + ' 


ib 


^l^m^n 


j(a;i)[j(a;m) • j(wn)].(13) 


where b = n^e'^b/jm^. 

It is easy to see that up to third order in the external 
electric field E the total current is given by: 

j(a;) = cr^^^(u;)E(w)-|-7o-^^^(w) ^ 

{Imn) 

(E(a;m) ■ E(a;„)), 
(14) 

where 7 = jm^jriQe^ and (T^^^(a;) is the linear conduc¬ 
tivity given by: 


1 -f 


ib 


^l^m^n 




uiaojuj) 
uj -I- if3ao{uj) ’ 


(15) 


Nonlinear optical wave generation 


The nonlinear conductivity leads to another interesting 
consequence when the optical response of the system is 
considered. We note that the spin liquid system can be 
treated as a nonlinear medium, where the linear electric 
susceptibility and the third-order nonlinear electric 

susceptibility are given by: 




iaoiu;) 


UJ + i/3(To(w) ’ 


(18a) 






1 -h 


ib 


^l^m^n _ 

(a;m)cr^^^ (w„).(18b) 


As shown in Ma and Ng [l^, the complex linear suscep¬ 
tibility leads to surface plasmon modes propagating along 
the interface between a linear medium and the spin liq¬ 
uids in the linear response regime. Assuming a boundary 
along y-direction at a: = 0 , the surface mode has a form 
E 5 (r, t) = -|- c.c., where ki = {kis,k,0) 

is the wave vector, and Ai(aj) is the amplitude of the 
electric field. The components of the wave vector ki are 
given by 0 : 


2 ^ ^_eo£77L 

Eo + ' 

2 ^ ^j£7m1!_ 

£0 -I- £(^)(w) ’ 


(19a) 

(19b) 


where cro(w) is the ac Drude conductivity cro(w) = 
CTo/(l — iojr) and /3 = gv"^. We note that there is no 
distinction between longitudinal and transverse currents 
in the limit q —0. The linear conductivity cr*^^)(a;) is 
consistent with the result obtained through loffe-Larkin 


where £0 is the dielectric constant of the linear medium 
(the vacuum), and £^^^( 0 ;) = 1 - 1 - 47 ry(^)(u;) is the dielec¬ 
tric function of the nonlinear medium (the spin liquids). 
We are interested at surface modes where k is nearly real 
and kis is nearly imaginary. It was shown in Ma and Ng 
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[l^ that a sizable frequency range where this condition 
is met is found in spin liquids (with large spinon Fermi 
surfaces) close to metal-insulator transition. 

The surface wave will induce a nonlinear electric po¬ 
larization: 

P3(3w) = (20) 


where 03 is the direction of the polarization vector, and 
k 3 s = 3ki. The nonlinear polarization P 3 ( 3 a;) acts as a 
source term in the nonlinear wave equation 


V X V X Es( 3 w) — 


{3uj) ^(i)( 3 ^( ^ Es(3a;) 
^^^%^e 3 P 3 ( 3 cc)e*'^^-^ ( 21 a) 


and generates radiation with frequency 3a; [20|, [2l| . The 
solution of the nonlinear wave equation contains a partic¬ 
ular part with wave vector kss and a homogeneous part 
with wave vector ks where fc| = (3aj), i.e., 

E«(3a;) = A 3 p{3u;) (3a;)e**^=*F (21b) 


In the presence of boundary, the particular and homo¬ 
geneous solutions are coupled. Solving the wave equa¬ 
tions with an incoming wave with frequency uj at angle 
9i upon the boundary at x = 0 (see Fig. |TL we find that 
the particular solution can be chosen as [ 20 , [ 2 l| : 


A3p(3a;) 


47rP3(3a;) / _ k3s(k3s ■ 63) \ 

es(3a;) — e(^l(3a;) V ^ )' 


where £s(3aj) = Matching the phase factors at 

the boundary leads to Snell’s law: 



FIG. 1. The surface modes (with frequency oj) supported by 
the linear complex susceptibility (Here the Otto configu¬ 
ration has been used to excite the surface modes: the three 
layers (1|2|3) are prism, linear medium (e.g., the vacuum) 
and the spin liquids, respectively.) will induce the waves with 
frequency 3a; that can propagate in the spin liquids. If the 
propagation angle 63 <C 7 r/ 2 , it can be distinguished from 
other surface waves and is detectable. 


For 03 <C 7 r/ 2 , the wave will propagate away from the 
boundary and can be distinguished from the other surface 
waves as shown schematically in Fig. |TJ 

For the 17(1) spin liquids, the linear dielectric function 
can be written as [l^ . 


k 3 ,y = kss^y = 3k. (22) 

and the amplitude of the homogenous part is found to 
be: 

^ _ 47rP3(3a;) y/so cos 03 s + \/£s (3a;) cos 0/ 

e,j(3a;) — e(il(3a;) cos 03 -I- y/£(^)(3a;) cos 0/ ’ 

where 03 ^ and 63 are the propagation directions of the 
particular part and the homogenous part, respectively. 
03 s is related to 0/ and uj through Snell’s law: 


e^^^(a;) = e(oo) 





2 

p 


+ iUJT 


-1 ’ 


(25) 


where £( 00 ) is the background dielectric constant, and 
tCg = / 3 a;p/ 47 r with ujp the plasma frequency. With 
Eq. (IMl) . it is straightforward to show that in the limit 
T—>- 00 , 03<7r/2if 


1 1 1 

£li)(3a;) £(i)(a;) ^ £0 ’ 


(26) 


sin 0 / 



(23) 


The particular part (with wavevector kss = 3ki = 
(3A:is, 3fc, 0)) decays rapidly in the nonlinear medium as 
kis is nearly pure imaginary when k is real. However the 
homogeneous part decays rather slowly and can prop¬ 
agate in the nonlinear medium [l^ . The propagation 
direction of the waves is given by 


tan 03 


k3,x 


3k 

^(^£(i)(3^) _ (3A:)2 


(24) 


where we choose £^^^(3a;) > 0. For nonzero r this condi¬ 
tion is satisfied in a narrower frequency region. The real 
and imaginary parts of 63 is shown in Fig. [2] for r = O.lw 
with a background dielectric constant £( 00 ) = 2.5. We 
see that there is a finite frequency range where Im 03 is 
small and Re 03 < 7 r/ 2 . In particular, the propagation 
angle of the waves in third harmonics generation with 
frequency 3a; is 0 ~ 52°, generated by an incident wave 
with frequency oj = 0.4a;p, which also lies in the fre¬ 
quency range where the surface plasmon mode exists. 

Lastly we estimate the magnitude of the nonlinear ef- 
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FIG. 2. Real and imaginary parts of the propagation angle 
63 of the waves with frequency Sw. We take t~^ = O.luj 
in the calculation [l8l| . It clearly shows that in some fre¬ 
quency region, the angle is nearly real 63 -C n/2. If we choose 
e(oo) = 2.5, for incident wave with frequency to = OAujp, the 
propagation angle is about 83 ~ 52°. 


Conclusion 

Summarizing, we study in this paper the nonlinear re¬ 
sponse of 17(1) spin liquids under the perturbation of a 
time-dependent but spatially uniform electric field within 
the framework of a 17(1) spin liquid theory. The third- 
order optical conductivity (w) is determined where we 
show that it has the same low energy power-law struc¬ 
ture as the linear optical conductivity (Tl^)(a;). We show 
that the nonlinear optical conductivity leads to an inter¬ 
esting, observable phenomenon where the surface modes 
with frequency ui supported by the complex linear suscep¬ 
tibility generates through the third-order process an EM 
wave with frequency 3a; which can propagate inside the 
nonlinear medium. This excited mode can be detected 
in an Otto configuration. 

We acknowledge support from HKRGC through grant 
No. 603913. 


feet. We consider the nonlinear susceptibility 


(3a;; a;, a;, a;) 


yX^^^So;) 



[cr(i)(a;)]^ 


In the typical frequency scale ojp ~ where the 

surface plasmon emerges, the linear conductivity lies 
in the scale ~ 0 “^cm“^ 18|, and the magni¬ 


tude of the linear susceptibility is of order x^^H^w) 
(3a;)/3a; ^ 10“^°fl“^m“^s. From Eq. (ITOl) . we 
also find 7 ^ l/(fc|.r) ~ r“^/( 2 mep), where ep ~ 
lOmeV is the Fermi energy. Choosing r ^ and 

tT^^)(a;)/(To 1 and 7 = = ye^r^/erg. The 

first term of the nonlinear susceptibility is thus of order: 




eV" 

2mep 


X^^^So;) 


(a;) 

CTo 


n 3 


10 " 


Vv^ 


For the second term, q is a dimensionless coupling con¬ 
stant of order unity 13| , and a simple dimension analysis 
gives 61 ~ Ijgv'lE'^, where Ec is a microscopic energy 
scale for charge fluctuation. For system close to metal- 
insulator transition, E^ ~ ep, and Va ^ vp ^ 10 “^c. 
Thus, as a rough estimation, the second term is of order 


X^^^(3a;) ~ b[x^^\uj)Y 






10 ' 


2 /y 2 ^ 


which is smaller than the first term in the frequency scale 
considered. We can see the magnitude of the nonlinear 
susceptibility is comparable with the typical value x^^^ ^ 
10“^°m^/V^ [2lj. Thus we expect the third harmonic 
generation process can be observed experimentally in the 
[/(I) spin liquids. 
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